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The general method for treating non-Gaussian wave functionals in the Hamiltonian formulation 
of a quantum field theory, which was previously developed and applied to Yang-Mills theory in 
Coulomb gauge, is generalized to full QCD. The Hamiltonian Dyson-Schwinger equations as 
well as the quark and gluon gap equations are derived and analysed in the IR and UV momentum 
regime. The back-reaction of the quarks on the gluon sector is investigated. 
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1. Introduction 

Coulomb gauge Yang-Mills theory has attracted a considerable amount of attention over the 
last few years. In the continuum, both the Hamiltonian ^, |3| [|] and the Lagrangian [Q, |6[ ^ ||] 
approach have been investigated. In the Hamiltonian approach, in particular, the use of variational 
methods with Gaussian wave functionals has led to the so-called gap equation for the inverse equal- 
time gluon propagator. The analytical and numerical solutions show an inverse gluon propagator 
which in the UV behaves like the photon energy but diverges in the IR, signalling confinement. 
The obtained propagator also compares favourably with the available lattice data [§]; deviations 
in the mid-momentum regime (and minor ones in the UV), which can be attributed to the gluon 
loop escaping the Gaussian wave functionals, can be taken into consideration by going beyond the 
purely Gaussian form. To treat these non-Gaussian functionals the authors developed in Ref. JTo| ] 
a method relying on Dyson-Schwinger Equations (DSEs) by exploiting the formal similarity be- 
tween vacuum expectation values in the Hamiltonian formalism and correlation functions in Eu- 
clidean quantum field theory. In this talk we report about the implementation of these techniques in 
the quark sector of the theory, in order to investigate issues like the spontaneous breaking of chiral 
symmetry. 

The topic of chiral symmetry breaking in Coulomb gauge has been studied for example in 
Ref. [p|]. While it has been shown that an infrared enhanced potential can account for chiral 
symmetry breaking, the calculated physical quantities, such as the dynamical mass and chiral con- 
densate, turn out to be far too small. It has been recently shown [ |12| ] that using a wave functional 
which includes the coupling of the quarks to the transverse gluon field improves the results to- 
wards the phenomenological findings. We apply here the techniques developed in Ref. []!(]] to the 



wave functional proposed in Ref. [12]. While we do not expect considerably different results from 



Ref. | |12| ] for physical quantities, the approach presented here has a broader range of applications, 



and could easily be applied to more complicated wave functionals. 



2. Hamiltonian Dyson-Schwinger Equations 

The derivation of the Hamiltonian DSEs in pure Yang-Mills theory has been presented in 
Ref. [0]. We summarize here briefly the essential steps leading from the choice of the wave 
functional to the pertinent DSEs. 

In the Schrodinger picture of the Yang-Mills sector we represent the gauge field and the canon- 
ical momentum operators in the basis |A) of the eigenstates of the field operator by 

(A\A\(j>)=A0[A}, ( A \n\<!>) = ^&[A], (2.1) 

where 0[A] = (A\<&) is a physical state. The vacuum expectation value (VEV) of an operator K is 
therefore given by 

(K[A,U\) = J @AJ A 4>*[A]K[A,^]<t>[A]. (2.2) 

In Eq. (2.2) the functional integration runs over transverse field configurations satisfying the Cou- 
lomb gauge condition, d,A" = 0, and is restricted to the first Gribov region; J a = Det(G^ 1 ) is the 
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Faddeev-Popov determinant which arises from the gauge fixing, and 



Km 



(2.3) 



is the inverse Faddeev-Popov operator, with g being the coupling constant and f ahc being the 
structure constants of the su(A/ c ) algebra. 

In a similar manner, a state \<P) in the fermionic Fock space possesses a "coordinate" repre- 
sentation 

{S\*) = *(S\Z), (2.4) 

where £, , E, ' are complex Grassmann fields. The VEV of a fermionic operator in this state is given 
by 

<^[^,^t^ = y n ^t^ e -^-(A + -A-)^*(^^t ) ^^^t] 0g\£) } (2.5) 

where the A± are the projectors onto states of positive and negative energy of the free Dirac theory. 
The Dirac field operators act on functionals according to 



8^ 



(2.6) 



We have put explicitly a hat over the fermion operators in Eq. {!£) to distinguish them from the 
Grassmann (classical) fields used in the "coordinate" representation. The exponential factor 
occurring in the fermionic functional integration in Eq. (1.5) arises from the completeness relation 
for fermionic coherent states, see e.g. Ref. [13]. 

The vacuum state of QCD can be assumed to be of the form 



1 



f[A, % , ^] =: exp <J - - S A [A] - S f [£, , A] 



(2.7) 



where Sa is a functional of the gauge field only, while Sf contains the fermion and fermion-gluon 
interaction parts. Dyson-Schwinger equations are derived by starting from the identity [10] 

= y'^A^t^^| /Ae -^(A + -A_)^* [A ^^t] jR: [ A ^^t]^[ A ^^t]j (28) 
with^€{A,^^}. 

To proceed further, we need an explicit ansatz for the vacuum wave functional Eq. (2.7). Since 
we are interested here mainly in the quark sector, we choose the simple Gaussian form 

S A = J d 3 x d 3 yA«(x) C0(x,y)Af(y) (2.9) 

for the Yang-Mills part. It should be stressed, however, that this is in no way necessary, and non- 
Gaussian generalization can be taken as well [10]; we restrict ourselves to the form Eq. (2.9) only 
for pedagogical reasons. 



For the quark wave functional we choose the ansatz [12] 



(2-10) 
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Figure 1: Hamiltonian DSEs for the quark (left) and gluon propagator (right). Filled dots stand for dressed 
propagators and empty circles for proper functions. The small empty boxes represent the variational kernels. 
Curly lines represent gluons, dashed lines ghosts, and straight lines fermions. 



The kernel contains both the purely fermionic part and the coupling of the quarks to the trans- 
verse gluons, 

Kr(x,y) = 8 mn p s(x,y) + g f d\V(x,y;z) a-A a {z) C (2.1 1) 

Here, j8 and a are the usual Dirac matrices, t a are the generators of the group in the fundamental 
representation (or, in general, in the representation to which the fermion fields belong), and the 
functions s and V are the variational kernels. 

With these choices for the wave functional the resulting DSEs are represented diagrammati- 
cally in Fig. [I[ 



3. Variational Approach to QCD 



The Hamilton operator of QCD in Coulomb gauge [|14J], resulting from the resolution of 
Gauss's law in the canonically quantized theory, can be written as 

+ J d^xf^x^l-iOidi + Pmlxir^-g J d\^ m '\x)aiA c f{x)C^(x) 

+ J d 3 xd 3 yJ A [ p a (x)J A Ff(x,y)p b (y). (3.1) 



In Eq. O 

F?j{x) = diA^x) - djA?(x) +gf' hc A l l(x)A){x) 
is the spatial part of the field strength tensor, 



Ff(x,y) = f d\G a A c (x,z) (-d?)Gf(z,y) 



(3.2) 



(3.3) 



is the so-called Coulomb kernel, which arises from the longitudinal component of the electric field, 
and 
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p a (x) = y m \x)t™ n xf{x) +f abc A b i {x 



■abc \ b i 



i8A\(x) 



(3.4) 



is the total colour charge density. 

The energy density can be evaluated as expectation value of the Hamiltonian Eq. (3J). The 
Hamiltonian DSEs in the general form Eq. (2^) enter this calculation twice: first to express the 
VEVs of the operators iff, in terms of the Grassmann fields t,, E, \ and second to turn the result- 
ing expressions, which involve propagators and vertex functions, in a functional of the variational 
kernels. In this calculation we have to introduce approximations, and we evaluate the energy den- 
sity up to two loops (see Ref. [10] for more details). 
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4. Gap Equations 

The variation of the energy density with respect to V fixes the vector kernel 

l+s(p)s(q) 



"(n I 7f\ I \ ^S 2 {P)+MPH4) i | ? | l-.v 2 (9)+2,(j?).v@ ' (4 - 1} 

U(p + q) + \p\ ITF p + |9| j^pgj 



where £2 = i(AA) 1 is the inverse gluon propagator. Replacing the kernels occurring on the right- 
hand side by their tree-level forms, i.e. Q.(p) = \p\ and s(p) = 0, we obtain 

V (p,q) = - ■ * , M , (4-2) 
\p + q\ + \p\ + \q\ 

which is exactly the leading-order perturbative expression for the quark-gluon vertex for massless 
fermions JT5|], 



Equation (4.1) can be inserted back into the expression for the energy density, and taking 
functional derivatives with respect to (o(p) and s(p) we obtain the gluon and quark gap equations. 
As a first approximation, we will ignore the ^-dependence of the denominator of Eq. (|Q1). Keeping 
the whole denominator structure would give rise to more complicated expressions which, however, 
have the same IR and UV asymptotic behaviour. 

The approximated gap equation for the scalar quark kernel s reads 

g 2 C F f d 3 q s(q)[l-s 2 (p)}-p.qs(p)[\-s 2 (q)] 
\P\ S \P) = ^ / 7KZ^F{p — q)- 



+ 



2 J (2?r) 3 u ^ l+s 2 {q) 
g 2 C F f d 3 q X(p,q) [l + s(p)s(q)] [s(q) - s(p)] 
2 J (2ny a(p + q) [l+s 2 (q)} [Cl(p + q) + \q\ + \p\] 



(4.3) 



where F = (Fa) is the Coulomb propagator and X is a tensor structure arising from the traces in 
Dirac space. The first term on the right-hand side of Eq. (Ob was obtained by Adler and Davis 



[11]; the second term is due to the coupling of the quarks to the transverse gluons. As mentioned 



in the Introduction, in Ref. [|12|] the coupling to transverse gluons was also taken into account. The 



authors of Ref. Q12[ ] used the same wave functional Eq. ( |2.10| ) as in this work but did a quenched 
calculation. As a consequence the effect of the Yang-Mills kinetic energy operator on the quarks 
ws neglected. 



Trading the scalar kernel s in favour of the mass function M, Eq. ( |4.3D takes the following 
simple form 

Cr t i\ „,_ _>|P|M(?)-MI?|M(P) 



+ 



(2k) 3 ^ ^ S(q) 
2 C F f d 3 q X(p,q) \p\M(q)-\q\M(p) (4 ' 4) 



2 J (271)3 Q(p+4)[a(p + $ + \p\ + \3\] S(q) 
with 



£(p) = Vp 2 + M 2 (p). (4.5) 



5 



Hamiltonian Dyson-Schwinger Equations of QCD 



Davide Campagnari 



0.5 
0.45 

0.4 
0.35 

0.3 

D(p) 0.25 
0.2 
0.15 
0.1 
0.05 






Gribov 




with quarks 


- f\ 

/' v\ 

- \\ 




I \\ 
- j \\ 




[ 

1 \\ 




1 N \ 
f 









3 4 5 
p/V*c 



Figure 2: Gluon propagator D = 1/ (2£1) in the presence of quarks for the gap equation ( 4.6 ) assuming the 
Gribov form for the Yang-Mills part. Physical dimensions are set by the Coulomb string tension a c . 



The variation of the energy density with respect to the gluon kernel ft) combined with the gluon 
propagator DSE yields the gap equation for inverse gluon propagator Q. 



p + Yang-Mills terms 

d 3 q [l+s(q)s(p + q) 



2Q.(p) + \q\ + \p + q\ 



(27T)3 [l +lS 2@)][l +lS 2(£ + 2)] [ n(/?) + 0| + |£ + 2|]2 

p-q+(p-q) 2 



l + - 



\q\ \p+q\ 



(4.6) 



where the usual Yang-Mills terms (ghost loop, Coulomb interaction, and possibly additional terms 
stemming from a non-Gaussian ansatz) can be found in Ref. 



Equations ( |4.3| ) and ( |4.6| ) form a coupled system, toghether with the DSE for the ghost propa- 
gator and the Coulomb form factor. The analysis of the full coupled set of equations is subject of 
ongoing work. As a first estimate, we investigate the effect of the quark loop on the gluon propa- 
gator: we take the Gribov formula for the pure Yang-Mills gluon propagator and use for the scalar 
kernel s a form fitted from the data in Ref. [12]. The result is plotted in Fig. ||. As we see, the 
gluon propagator loses some stregth in the mid-momentum regime. This behaviour is known from 
Landau gauge studies, and Coulomb gauge investigations on the lattice confirm this JlB]]. 

5. Conclusions 

The general method to treat non-Gaussian wave functionals in the Hamiltonian formulation of 



a quantum field theory presented in Ref. [JlOp has been applied here to QCD in Coulomb gauge. 



We have used the quark wave functional suggested in Ref. []12[], which includes the coupling of the 
quarks to the transverse gluons. The resulting quark gap equations are similar in both approaches; 



the numerical solution of the equations and the comparison to the findings of Ref. [12] are the 
subject of ongoing work. As a first application we have presented here the effect of the back- 
coupling of the quarks on the gluon propagator. In future work we intend to expand the present 



approach to QCD at finite temperature [17] and baryon density. 
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